The dynamics of azurite Cus (003)2(011)2 in a magnetic field as determined by 

neutron scattering 
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Azurite, a natural mineral made up of CuO chains, is also an intriguing spin- 1/2 quantum mag- 
net. There has been much debate as to the 1-dimensional (ID) nature of this material by theorists 
and experimentalists alike. The focus of this debate lies in the interactions between Cu-ions within 
the antiferromagnetically ordered state below 1.9 K. We present high-resolution inelastic neutron 
scattering data which highlight the complexity of the magnetic ground state of azurite. The appli- 
cation of magnetic fields and temperatures were used to probe the excitations revealing important 
information about the dynamics in this system. From this we are able to conclude that the ID 
Heisenberg antiferromagnetic spin chain model is not sufficient to describe the dynamics in azurite. 
Instead additional coupling including interchain interactions and an anisotropic staggered field are 
necessary to fully model the observed excitations. 

PACS numbers: 75.25.-j 75.50.Ee 61.05.fm 75.30.Et 



I. INTRODUCTION 



Low-dimensional magnetic systems, in particular cop- 
per oxides, are interesting subjects of study due to the 
novel physics that can emerge when they are cooled to 
low temperatures. Magnetism in these materials arises 
from the spin-1/2 moments carried by the Cu^+ sites, 
and the magnetic couplings can vary strongly in mag- 
nitude and sign depending on the geometry and en- 
vironment of the copper oxide bonds involved. Azu- 
rite, Cu3(C03)2(OH)2, is one such copper oxide material 
which has been the subject of debate in recent years 
A natural mineral, azurite can be considered one of the 
first experimental realizations of the ID distorted dia- 
mond chain model. While it orders antiferromagnetically 
below 1.9K^'^^ there has been much debate involving the 
relative exchange interactions between the Cu ions as 
shown in Fig. [TJ Much experimental analysis and theo- 
retical modeling have been performed however there is as 
yet no general consensus as to the exact magnitudes of 
the exchange interactions either within the ID diamond 
chain, or between chainsi"— i^. Despite this, there is sig- 
nificant evidence that suggests the exchange interaction 
J2 is by far the strongest interaction leading to a cou- 
pling of Cu-ions into dimers and monomer ^i^iii. From 
these results it is now widely believed that azurite can- 
not be described in terms of a simple isotropic exchange 
chain Hamiltonian, but rather that inter-chain coupling 
and/or anisotropic exchange must also be taken into ac- 



count when describing this material. 

Recently detailed investigations of the low temperature 
nuclear and magnetic structures has been performed ^^-s^. 
From these studies it was found that the symmetry of 
azurite was P2i rather than the previously determined 
P2i/c. For the first time, low temperature lattice pa- 
rameters were found at 1.28 K to be a = 4.99995(11) 
A, 6 = 5.82256(14) A, c = 10.33723(19) A with p = 
92.2103(17)° The magnetic moment structure was also 
revealed in this work where two different moments were 
found: mo = 0.684/i^ on the monomer site (Cul) and 
mo = 0.264/i^ on the dimer site (Cu2). From this work 
one can now speculate how the relative magnetic moment 
orientation of the dimer and monomer spins can affect 
the local magnetic environment in the ordered phase. 
For instance, it was proposed that the non- negligible 
moment on the dimer sites originates from the effective 
staggered field from the antiferromagnetically coupled 
monomers.—. 

In this paper, we present inelastic neutron scattering 
data of the quasi-one-dimensional (ID), S = 1/2 anti- 
ferromagnet (AFM) Cu3(C03)2(OH)2. These investiga- 
tions focus on the magnetic excitations in azurite for 
applied fields below the plateau phase - that is below 
IIT when applied perpendicular to the Cu-chain direc- 
tion. We will compare our data with the ideal Heisenberg 
model and reveal the inconsistencies between the two. 
We use the recent structural and magnetic information 
to devise a model which better describes azurite. 
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FIG. 1: (Colour online) Diamond chain model showing the 
structure of azurite and the relative exchange interactions 
within the chain. It is now widely believed that J2 is the 
strongest interaction in azurite. 



II. EXPERIMENTAL METHOD 

Inelastic neutron scattering measurements were con- 
ducted at the Helmholtz Zentrum Berlin (HZB), us- 
ing FLEX, the cold-neutron triple-axis spectrometer. 
A large, naturally obtained single crystal sample of 
Cu3 (003)2 (OH) 2 which has been studied previously^ was 
used for all inelastic neutron scattering measurements. 
The sample was aligned with either an a6-scattering 
plane or a 6c-scattering plane such that measurements 
could be taken along each of the principal lattice di- 
rections i7, K and L. The sample was kept at 1.5 
K to ensure that the system was in the ordered mag- 
netic phase. Oonstant-Q, energy-transfer scans were per- 
formed within the range 0-2.5 meV with incident neu- 
trons fixed at ki = 1.3 A~^. The collimation was set to 
guide-60 '-open-open with a horizontally curved analyser. 
This combination of parameters gave an instrumental en- 
ergy resolution of 0.11 meV. 

The high resolution time of flight spectrometer OSIRIS 
at ISIS, UK was also used to probe the magnetic exci- 
tations in azurite. Magnetic fields up to 7T were ap- 
plied perpendicular to the chain direction (i.e. perpen- 
dicular to the crystallographic 6* axis). The sample was 
mounted onto the end of a Kelvinox dilution insert and 
placed within the 7.5 T vertical field magnet. The tem- 
perature of the sample was varied between 0.05 - 20 K 
but was stable for each measurement to within ±0.01 K. 
The sample alignment was checked in diffraction mode, 
using the backscattering detector banks confirming an a* 
- 6* scattering plane. Data were taken using the chopper 
package PG002 which covered an energy transfer range 
of —0.2 < E < 2.5 meV with a resolution at E = of 
0.025 meV. This energy range was chosen to maximize 



the resolution of the AFM spin chain which was previ- 
ously observed below 2.5 meV energy transfer^. 



III. RESULTS 

A. Dimensionality of the dynamics 

Constant- Q energy scans were taken along each of the 
principal lattice directions K and L to ascertain the 
ID nature of the low energy excitations within the or- 
dered magnetic phase. These measurements were taken 
along the [H, 0.5, 0], [1, K, 0] and [0, 0.75, L] direc- 
tions and can be seen in Fig. [2l Since recent theoretical 
work has implied significant interchain contributions to 
the interactions in azurite, it was important to measure 
the dispersion relation along each of the three directions. 
From the data in Fig. [21 there is little or no observed 
Q-dependence along the H and L directions, within the 
resolution of these measurements. In contrast to this, 
there is a clear Q-dependence for measurements taken 
along the i^-direction - that is parallel with the diamond 
chains of Cu-ions. This is in fact the same result stud- 
ied and has been attributed to the coupling of the 
monomer ions along the diamond chain. These results 
imply that, within the resolution of these measurements, 
the interchain couplings along the L direction are at least 
an order of magnitude less than along K. 




0.25 0.5 0.75 
(0, 75. L] reciprocal lattice units 



0.00 0-25 0.50 0.75 
[0,0.75.L] reciproGal lattice units 



FIG. 2: Constant- Q energy scans taken along each of the 
principal lattice directions i7, K and L in azurite (labeled 
A, B and C respectively). The intensity of each data set 
was normalized to the elastic intensity. Black spots indicate 
the point density while smoothing has been applied. The 
dashed line in B is a guide to the eye. D shows a cosine fit to 
the excitation observed in C revealing an amplitude of 0.03 
meV. This indicates minimal dispersion along the L-direction 
around an order of magnitude less than along the i^-direction. 



3 



B. High resolution TOF 

The dynamics within the HKO plane were measured 
on the single crystal sample by summing together data 
from 20 different sample rotation positions covering 90°. 
These results reveal the rich nature of the low energy 
magnetic scattering excitations in azurite. It is believed 
that this scattering results from an effective ID chain 
primarily of monomer spins, while the dimerised pairs 
of Cu-ions primarily contribute to scattering at higher 
energies. In this way multiple Brillouin zones could be 
mapped out showing the dispersive relation along the 
chain direction, as shown in Fig. [3l These data have been 
treated for detector efficiency with a background sub- 
tracted and the intensity modulation across the Brillouin 
zone is a result of the Cu-form factor. For these measure- 
ments, the sample temperature remained constant at 60 
mK in order to reduce the influence of thermal fluctu- 
ations on the excitations. While the overall dispersion 
looks very similar to that observed on FLEX and in ear- 
lier measurements^, there are clearly additional features 
which are revealed from the higher resolution study. In 
particular, there is evidence now of 2 gaps at the AFM 
Brillouin zone center (0, -0.5, 0) located at Ai ^ 0.4 
meV and A2 ~ 0.6 meV. While the periodic dispersion 
resembles the scattering from a ID Heisenberg Antifer- 
romagnetic chain (HAFC) with a spinon continuum, this 
model is gapless as demonstrated in the simulation of 
Fig. 4c. Heisenberg interchain interactions would al- 
low for the observed gaps in the iC-dispersion however 
these would also be responsible for dispersive behaviour 
along the other principal directions, H and L. Since Fig. 
[2] showed almost no dispersion, then we must consider 
an anisotropic interchain interaction which would not be 
dispersive perpendicular to the chain. 




-2 -1,5 -1 -0.5 
[0, K, 0] reciprocal lattice units 



C. Temperature variation of the low energy band 

Focusing on one Brillouin zone with higher statistics, 
the temperature dependence of the low energy scatter- 
ing can be observed, as seen in Fig. HI The two gaps 
at the zone center are clearly evident in the lowest tem- 
perature data set, however at higher temperatures they 
are not observed. At 2.5K azurite is no longer in its 3D 
magnetically ordered phase. However scattering from the 
diamond chain and the spinon continuum is still evident. 
In particular the scattering is strongest at the zone center 
with a possible sinusoidal dependence resembling that of 
the IDHAFC model. It is clear that while the signal is 
much weaker at 2.5K, the chain dynamics still persist. 
While the counting statistics of the higher temperature 
data sets are insufficient for a quantitative assessment, 
they show a qualitative agreement with the ID spin chain 
model. 




-1 -D.8 -0.6 -D.4 -D.2 D -1 -0.6 -O.S -0.4 -0.2 

[0, K, □] reciprocal lattice unrt [0. K,OJ reciprocal Mice units 




|OjK,0| reciprocal lattice units [Ci,K,0] reciprocal lattice ufifts 



FIG. 4: The dynamic scattering from azurite along the Cu- 
chain direction in azurite is shown at 0.04 K (A) and 2.5 K 
(B). The periodic dispersion in K at 40 mK resembles the 
ID HAFC with a spinon continuum however this model does 
not allow for the energy gaps, Ai and A2, clearly observed 
in this high resolution data (indicated by the arrows). The 
lower plots show the expected scattering from a ID HAFC 
at 0.04 K (C) and 2.5 K (D). The absence of gaps from the 
model in C indicates that interchain coupling and anisotropy 
terms play a significant role in the azurite Hamiltonian. In 
contrast the data and simulation taken above the magnetic 
ordering temperature look quite similar. 



FIG. 3: Inelastic neutron scattering data covering two Bril- 
louin zones showing the periodic dispersion along K. Mea- 
surements were taken at 60 mK. 



D. Applied field TOF 



By applying a magnetic field we can influence the mag- 
netic dynamics in such a way as to reveal more about the 
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interactions within azurite. The dynamics of azurite were 
observed as a function of magnetic field up to 7T apphed 
parallel to the c-direction. The zone center and bound- 
aries were observed to show the most notable changes 
and have been plotted as a function of field at 0.04K in 
Fig. [SI The zone-boundary mode at a zero field energy 
of 1.2meV clearly splits into 2 modes at low fields with a 
third mode splitting from around 3T. While these modes 
diverge with increasing field, they do not follow a Zeeman 
type splitting. 

At the zone-center, three modes are clearly visible with 
varied behaviour in field. The weak upper mode at ^ 1.2 
me V tracks the zone boundary behaviour of Fig [5] A and 
could be a shadow due to secondary scattering processes 
from incoherent scattering. The lower gap, Ai, shows a 
clear decrease in energy with applied field while the upper 
gap, A2, remains approximately constant with field. The 
different field behaviors indicate that the origin of the two 
gaps also differs as will be discussed later. 
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FIG. 5: The field dependence of the low energy excitations 
taken at the zone boundary (A) with Q = [0, —0.25, 0] and 
the zone centre (B) with Q = [0, —0.5, 0]. Gaussian fits to the 
excitations show a divergence of the modes at the zone center 
(C) while each of the three modes at the zone boundary (D) 
show different field dependent behaviour. These data indi- 
cate that the nature of the excitations is not that of standard 
spinwaves in an isotropic antiferromagnet. 

The effect of the applied field on the Q-vector of the 
energy gaps can reveal information about the character 
of the two different modes. Fig. [6] shows how the lat- 
tice vector of the dynamics behaves in fields up to 7 T. 
The lower gap, which has been displayed by binning the 
energy range 0.35 < En < 0.45 meV, shows almost no 
deviation from the zero field value. The gap at higher 
energies 0.55 < En < 0.65 meV shows a different field 
dependence with a clear splitting of K from around 1 T. 
As discussed below, this suggests that the lower mode 
has originated more from the transverse continuum (spin 
fiuctuations perpendicular to the field) and the upper 
from the longitudinal continuum (spin fiuctuations along 
the field). 



0.35 < E < 0.45 meV ^ 


^^KE<a^5 meV 
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FIG. 6: The variation in ii^- vector of the zone-center energy 
gaps with applied field. Left: the lower gap located in the 
range 0.35 < En < 0.45 meV shows almost no change in 
K with field (as indicated by the black line) until 6T when 
Gaussian fits indicate a small splitting. Right: the upper gap 
in the range 0.55 < En < 0.65 meV shows a more significant 
splitting at low fields. The black lines follow the magnetiza- 
tion per spin of the Heisenberg chain as defined by equation 

(1). 



IV. ANALYSIS AND DISCUSSION 

A. Comparison to the ID HAFC model in a field 

To gain more insight into the physics of azurite we 
first consider the effect of applied field and how this in- 
fiuences the spin 1/2 Heisenberg chain model. In Fig. [71 
A, C, and E we demonstrate how applied field will af- 
fect the spin correlations for the ID Heisenberg model. 
These simulations involve taking the average contribu- 
tion from the spin correlations along the field direction 
S^^{Q^uj) and perpendicular S^{QjUj). The spin correla- 
tions are computed for finite chains of 192 sites using the 
integrability-based ABACUS method^"i^. Overall, the 
general features observed in S{Q,uj) compare favorably 
with the neutron scattering data shown in Fig. [3B, D, 
and F despite the absence of gaps. However it is clear 
that the rotational symmetry in azurite is broken by the 
applied field. 

In zero field there is only one contribution to the scat- 
tering S{Q,uj) for the Heisenberg model which is shown 
in Fig. [8] A. However in an applied field, the separate con- 
tributions from S^^{Q^uj) and 6'^(Q,cj) combine to give 
the scattering pattern in Fig. [8] B for 4.5 T. Here, we 
define S^^{Q^uj) as the spin correlations longitudinal to 
the applied field and S-^{Q^uj) as the average of the spin 
correlations transverse to the applied field. In contrast 
to zero field, the 4.5 T scattering from the transverse and 
longitudinal correlations can be considered separately as 
in Fig. [8] C and D respectively. 

The transverse and longitudinal continua behave very 
differently in applied fields and have been modeled in 
Fig. [9l Here the transverse (left panels) and longitudi- 
nal (right panels) continua have been calculated for an 
anisotropic Heisenberg Hamiltonian, as defined in Equa- 
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0.6 0,8 1 -0.75 -0.5 -0.25 
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FIG. 7: A, C and E show simulations of the inelastic neutron 
scattering data for applied fields 2T, 4.5T and 6T. These can 
be compared to the experimental results for the same fields 
in B, D and F. A background from the zero field data set has 
been subtracted from the data on the right. 




0.4 0,6 O.a 1 0.2 0.4 0.6 
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FIG. 8: Simulation of the ID HAFC at zero field (A) and at 
4.5T (B). The 4.5 T simulation can be separated into the two 
components: transverse (C), indicating the spin correlations 
perpendicular S'^{Q,(jj), and longitudinal (D), indicating the 
spin correlations along the field direction S^^{Q,uj). 



tion 1 from Ref.^^, with applied fields increasing up to 
Hz = 0.1. Here the saturation field along (equiv- 
alent to 11 T for azurite) is Jx^ Jy and Jz dependent. 
The anisotropy in this XYZ model is responsible of the 
energy gap observed at the zone centre. As can be seen, 
the transverse mode does not split in K at the zone centre 
whereas the longitudinal does. However both show sig- 
nificant splitting at the zone boundary. Since the lower 
mode observed in azurite (0.4meV) (Fig. [6j) does not 
appear to split in K with field we can deduce that this 



originates from the transverse continuum while the up- 
per gapped mode (0.6meV), which does split with applied 
field, originates from the longitudinal continuum. 
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FIG. 9: Simulation of the anisotropic Heisenberg Hamilto- 
nian, as defined in Equation 1 from Ref.— , where = 1, 
Jy = 0.5 and = 1. The transverse (left column) and 
longitudinal (right column) continua scattering have been 
plotted for applied magnetic fields (from top to bottom) 
Hz = 0.01, 0.05, 0.1. While the zone centre of the transverse 
mode shows no incommensuration along K, the longitudinal 
mode shows significant splitting. Both show a strong field 
dependence at the zone boundary. 



The magnetization per spin in a spin- 1/2 Heisenberg 
chain versus field can be approximated using the follow- 
ing equation from Ref.^-: 



a(h) = —a sin ■ ; — - 

^ TT l-7T/2^7Thsat/2h 



(1) 



where < h < hsat and hgat = 11 T for azurite. 
The incommensuration along K (to be compared to the 
right panels of Fig. [9] for the longitudinal mode) is 
q(h) = 1/2 ± (7{h). This uses the relationship between 
magnetization and filling factor for spinous. It also in- 
cludes the field dependent quantum renormalization fac- 
tors. The incommensuration should go from 0.5 in zero 
field to 0.35 and 0.65 by 7 Tesla in azurite and be ap- 
proximately linear in this field region. This looks similar 
to what is actually observed in azurite in Fig. [6] and 
confirms that the approximation of Eq. 1 is sufficiently 
accurate with respect to the exact solution as seen in Fig. 

m 

To see more clearly how the continuum boundaries 
shift with applied field, a simulation of the inelastic neu- 
tron scattering data is plotted in Fig. [11] showing the 
field dependence of all modes together and taken as a 
cut along K between 0.22 and 0.28. This calculation 



6 



\0 0,25 




FIG. 10: (color online) Plot of the magnetization curve for 
the isotropic XXX antiferromagnet. The blue line is the ex- 
act curve obtained from the Bet he Ansatz. The green dashed 
curve is the approximation in Eq. 1. The subplot gives the 
difference between exact and approximate curves, which re- 
mains small throughout the parameter range considered. 



captures important features of the modes that were seen 
in Fig. [5l From the mode splitting at the zone boundary 
one can deduce that the character of the observed exci- 
tations is not that of standard spinwaves in an isotropic 
antiferromagnet. In spinwave theory the excitations are 
oscillatory modes which are not changed significantly at 
the zone boundary. The strong zone boundary behaviour 
with field is often seen for spinon systems with strong 
quantum fluctuation^. 
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FIG. 11: Simulation of the inelastic neutron scattering data 
from Fig. [5] A at the zone boundary in applied fields up to 
7T. 

It is clear from the current high-resolution results that 
the simple ID Heisenberg antiferromagnetic chain model 
is insufficient to fully describe azurite. When determin- 
ing a suitable model to describe the excitation spectra, 
additional terms such as anisotropy and interchain cou- 
pling must be considered. Despite the broad similarities 
between the ID HAFC model and the azurite data, these 



additional terms are more effective at describing the de- 
tails in the neutron scattering data. 



B. Anisotropy 

As observed in Fig. [9l by introducing anisotropy into 
the system, an energy gap opens up in the excitation 
spectra. This can be modeled by the anisotropic XYZ- 
Hamiltonian which can be written as: 



.ii^^islst,, + s-s-,,) (2) 



where and Sf_^i are the spin matrices at neighbour- 
ing sites and the exchange interactions Jx 7^ Jy ¥^ Jz- 
Here, the excitations in the pure XYZ chain are gapped 
spinous which form an empty particle and filled hole 
band^^. Neutrons measure continua of pairs of these 
excitations just as in the isotropic case. Although the 
spherical rotational symmetry is completely broken there 
is only one type of spinon which is now gapped. Within 
a mean field theory the particle and hole dispersions are 
given bySi: 



uj{k) = ±J{tcos (k) - hf + A2 sin^ ik) 



(3) 



where t and A are calculable from the Jr 



Jy and Jz 



and h is the magnetic field. The transverse continuum of 
scattering is given by pair-states of momenta ki and /c2 
at g = /ci + /c2+7r and energy uo = uj{ki)-\-uj{k2). Whereas 
the longitudinal continuum moment is q = ki — k2 with 
the same energy relation uj = uj{ki) -\-uj{k2). The energy 
gap for zero field is the same for the longitudinal and 
transverse continuum so does not explain the different 
gaps in the continua that we actually observed. 



C. Interchain coupling 

While anisotropy is likely to be responsible for open- 
ing one energy gap in the ground state excitation spec- 
trum, interchain coupling can be used to describe the 
multiple gaps observed. The effects of interchain cou- 
pling in azurite have previously been down played due to 
the minimal dispersion observed in the earlier inelastic 
neutron scattering results^, however more recently the 
importance of interchain coupling has been highlighted 
as necessary to fully describe the interactions between 
Cu-ions.^'^ In fact, some form of anisotropic interchain 
coupling may also account for the slight dispersion be- 
tween the chains, observed in Fig. [2l For KCuFa, a 
model HAFC with weak interchain coupling, the quan- 
tum field theory works almost perfectly and is accurate to 
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a few percent. In contrast, for azurite the agreement with 
the quantum field theory for coupled Heisenberg chains 
is only within 30%. This means that an anisotropic in- 
terchain coupling (Ising like) is a credible possibility and 
interchain couplings should be considered to have a sim- 
ilar infiuence to within-chain anisotropy. 

Interchain coupling theories which have been devel- 
oped for the Heisenberg case, may also be applied for 
the case of Heisenberg chains coupled through an Ising 
interchain interaction. These theories can give a rea- 
sonable approximation of interchain coupling involved 
in azurite. The formula can be adapted to explain the 
Neel ordering temperature Tn = 1.9 K, ordered moment 
{S) = 0.684/iB and gap energy A = 0.40 meV and can 
all be related to the interchain coupling J_\_ by using 
quantum field theory. For instance, we can estimate the 
interchain coupling by considering the relationship be- 
tween the energy gap and the ordered monomer moment 
as outlined in Essler et al.,^^. A combination of equa- 
tions 15 and 56 from this paper gives the simple relation 
|Jx| = l/J(0.175[A/mo])^ = 0.048 meV. Similarly, tak- 
ing the approach of Schulz^^, the interchain interaction 
can be approximated as 0.065 meV for a gap A = 0.4 
meV. If we consider the critical temperature, Tn = 1.9 
K, to infiuence the strength of the interchain coupling, 
we find that for a cubic lattice the interchain coupling 
would be \J±\ = 0.083 meV as outlined in the work of 
Bocqueti^^. A more sophisticated analysis by Irkhin and 
Katanin^^ gives a slightly improved estimate of the inter- 
chain coupling as \J±\ = 0.086 meV with the following 
formula mo ^ 1.085(7^/7)^/^ 

The effect of three dimensional ordering of an array 
of spin chains is to induce an effective staggered field on 
the chains. In azurite, the spinous in the continuum are 
bound into well defined modes by the interchain coupling. 
Within this spinon picture for the Heisenberg chain, the 
staggered field produces two modes, one with gap Ai 
which is a doublet which has a polarization perpendic- 
ular (transverse) to the staggered field (ordering direc- 
tion). A breather mode (consisting of 2 bound solitons, 
or 4 bound spinous) also forms with a gap A2 which has a 
polarization along the staggered field direction (longitu- 
dinal). In the special case of pure Heisenberg interaction 
the gap ratio Ai:A2 is l:^/3. In general this ratio will 
vary with anisotropy for an XYZ Hamiltonian. For the 
case where the applied field and staggered field directions 
coincide then the lower mode with gap Ai which has a 
transverse polarization will remain commensurate as is 
observed, whilst the upper mode with A2 will split in Q 
with field as also observed. 

In the current experiments the applied field was along 
the c-axis when the magnetic ordering was observed for 
the monomers to be in the ac plane at an angle of 55° 
from (M. The spin fiuctuations along the 6-axis then are 
both transverse to the field and to the magnetic order. 
They then should not split in q and will have a lower 
gap energy that changes only weakly with field. Neutron 
measurements are sensitive to the fiuctuations perpendic- 



ular to Q and since the chain direction {b) was aligned 
almost completely perpendicular to the wavevector trans- 
fer, the 6-axis mode (transverse) is more strongly probed 
rather than the a correlation. The correlations for the 
c component (longitudinal) will always appear strongly 
because they are vertical and therefore always perpendic- 
ular to q transfer. These will pick up the incommensura- 
tion from the longitudinal nature of the field in relation 
to the ordering and so we expect incommensuration of 
the upper mode as is seen. 

In the spin chain material CDC, which is a funda- 
mentally Heisenberg system with very small interchain 
coupling^^, similar behaviour to azurite is seen. In par- 
ticular solitons and breathers are observed as well as the 
incommensurate behaviour of these modes. However the 
situation is not identical to azurite which is probably an 
XYZ chain with larger interchain coupling. Also, in CDC 
the staggered field is due to the applied field combining 
with the Dzyaloshinskii-Moriya interaction whereas for 
azurite the staggered field is probably from an interchain 
coupling effect. 

Finally, the most powerful method to determine the ex- 
changes in a quantum magnet is to use the method of sat- 
uration (which was developed to understand CS2CUCI4 
ir^) . So far only the exchange couplings along the chain 
direction in azurite have been observed by applying fields 
above 11 Tesla^. It is interesting to note that the energy 
gap at the zone center does not close at the onset of 
the plateau phase {Hd ~ 11 T). This is a signature of 
anisotropy and/or interchain coupling. The scale of this 
gap at 11 T was around 0.142 meV giving an interchain 
coupling of J± ~ 0.071 meV which agrees well with the 
values deduced above. In addition the scattering at the 
zone center, q = (1,0,0), when linearly extrapolated as 
a function of field (as seen in Fig. 4 fron>^) also does 
not close at T. A possible explanation for this involves 
interchain coupling from the (1,0,0) that is ferromagnetic 
along 110, 100, and 1-10 directions. This requires more 
detailed analysis but may indicate the presence of frus- 
trated interchain interactions. To understand this sat- 
uration phase better, the dispersions along the a and c 
directions should be measured such that the couplings 
within the plateau phase can be extracted. 



V. CONCLUSION 

In summary, high resolution inelastic neutron scat- 
tering results have shown that azurite is not well de- 
scribed by the ID HAFC model. Instead, more com- 
plex interaction pathways are a necessary addition to the 
spin-chain Hamiltonian to describe the observed features. 
Anisotropic interchain interactions of around 10% of the 
intrachain interactions are required to explain the two 
gapped modes in the dispersion along the chain direc- 
tion. The 3D long range magnetic ordered state can also 
be attributed to the presence of these interchain interac- 
tions. 
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